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ABSTRACT  One-dimensional  (ID)  zigzag  [b^bl-phenyl-Cgrbuty- 
ric  acid  methyl  ester  (PCBM)  nanoribbons  are  produced  by  folding 
two-dimensional  ultrathin  PCBM  nanosheets  in  a  simple  solvent 
process.  The  unique  ID  PCBM  nanostructures  exhibit  uniform  width 
of  3.8  =b  0.3  nm,  equivalent  to  four  PCBM  molecules,  and  lengths  of 
20-400  nm.  These  nanoribbons  show  well-defined  crystalline 
structure,  comprising  PCBM  molecules  in  a  hexagonal  arrangement 
without  trapped  solvent  molecules.  First-principle  calculations  and 
detailed  experimental  characterization  provide  an  insight  into  the 

structure  and  formation  mechanism  of  the  ID  PCBM  nanoribbons.  Given  their  dimensions  and  physical  properties,  we  foresee  that  these  nanostructures 
should  be  ideal  as  acceptor  material  in  organic  solar  cells. 
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Fullerene  Ceo  and  its  derivatives  are 
extensively  used  in  organic  electronic 
and  photovoltaic  (OPV)  devices.^ 
One  example  of  such  a  derivative  is  the 
[6,6]-phenyl-C6i -butyric  acid  methyl  ester 
(PCBM),  a  Ceo  fullerene  with  a  chemically 
bonded  functional  group.  The  addition  of 
the  functional  group  increases  the  solubility 
of  the  Ceo  in  polar  solvents,  which  facilitates 
device  fabrication  via  a  solution-based  pro¬ 
cess.  The  functional  group,  on  the  other  hand, 
decreases  the  fullerene  symmetry  and  conse¬ 
quently  affects  its  crystallization.^  Although 
growth  of  crystalline  Ceo  micro-  and  nano¬ 
structures  is  frequently  reported,^“^^  only  few 
studies  describe  the  synthesis  of  crystalline 
PCBM  nanostructures.^'^  This  is  explained 
by  the  lower  symmetry  of  PCBM  molecules 
and  their  flexible  functional  groups,  which 
promote  a  more  complex  crystal  packing 
and  inhibit  the  possibility  to  tune  the  grown 
structures  to  different  morphologies.® 
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One-dimensional  fullerene  (Ceo)  struc¬ 
tures,  namely,  nanorods  and  nanoribbons, 
are  of  particular  interest  for  applications  in 
integrated  electronic  circuits  and  organic 
photovoltaic  devices  because  they  can  po¬ 
tentially  be  used  as  transport  channels  with 
high  electron  mobility^'^^  and  concurrently 
provide  high  contact  surface  area  versus  the 
other  constituent  in  the  bulk  heterojunction 
OPV.  However,  practical  examples  of  one¬ 
dimensional  fullerene  structures  are  limited 
because  common  crystallization  techniques, 
such  as  slow  evaporation,^ templating,^® 
vapor— solid  processes,^ ^  fast  solvent  evapo¬ 
ration,^®  and  liquid-liquid  interfacial  preci¬ 
pitation  (LLIP),^^'^^“^®  provide  nanorods  with 
average  diameters  larger  than  150  nm.^^This 
is  much  larger  than  a  typical  thickness  of 
the  OPV  film,  and  since  ideal  ID  structures 
should  be  of  the  same  magnitude  as  the 
exciton  diffusion  length  of  the  charge  car¬ 
riers  in  OPVs  (/.e.,  10  nm),^^  there  is  a  large 

VOL.  XXX  ■  NO.  XX  ■ 


*  Address  correspondence  to 
thomas.wagberg@physics.umu.se. 

Received  for  review  August  10,  2015 
and  accepted  September  18,  2015. 

Published  online 

1 0.1 02 1  /acsnano.5b04972 


©XXXX  American  Chemical  Society 

000-000  ■  XXXX 

www.acsnano.org 


Figure  1.  TEM  micrographs  of  (a)  thin  PCBM  nanosheets  produced  by  the  LLIP  method,  where  the  interface  is  disturbed  by 
ultrasonication,  (b)  PCBM  nanosheet  during  folding,  and  (c)  lower  magnification  showing  various  PCBM  scroll-like  structures. 


motivation  to  controllably  grow  suitable  nanostruc¬ 
tures  in  this  size  regime. 

Only  recently,  Maeyoshi  etal?^  reported  the  forma¬ 
tion  of  Ceo  (and  different  Ceo  derivatives)  nanowires, 
with  an  average  diameter  of  8  nm,  via  high-energy  ion- 
beam-induced  Ceo  polymerization,  and  indeed,  such 
structures  showed  good  performance  in  OPV  devices. 

Transformation  of  2D  systems  into  ID  nanostruc¬ 
tures  has  already  been  shown,  for  example,  by  scrolling 
graphene  sheets,^^“^^  and  inspired  by  such  ap¬ 
proaches,  we  here  describe  a  facile  method  to  synthe¬ 
size  highly  homogeneous  ID  PCBM  nanostructures 
(width  ~4  nm)  by  folding  thin  PCBM  nanosheets.  Prior 
to  the  folding  step,  PCBM  nanosheets  are  produced 
by  LLIP  method.  The  method  provides  a  low-cost  and 
scalable  technique  to  synthesize  small  PCBM  nano¬ 
structures  for  practical  application  in  organic  electro¬ 
nics  and  OPV  devices.  By  extensive  experimental 
characterization  and  theoretical  modeling,  we  exam¬ 
ine  the  structure,  morphology,  and  the  energetic  sta¬ 
bility  of  the  synthesized  structure  and  show  that  the 
folding  of  the  nanosheets  leads  to  the  formation  of 
homogeneous  PCBM  nanoribbons  with  well-defined 
zigzag  edges  and  that  these  nanoribbons  subse¬ 
quently  can  stack  together  in  bundles,  forming  nano¬ 
rod  structures  with  diameters  of  10—20  nm. 

RESULTS  AND  DISCUSSION 

Ultrathin  nanosheets  are  synthesized  by  the  LLIP 
method  at  the  interface  between  a  solution  of  PCBM  in 
chloroform  and  methanol,  where  the  interface  is 
strongly  disturbed  by  ultrasonication.  Figure  la  shows 
a  transmission  electron  microscope  (TEM)  image  of 
the  synthesized  PCBM  nanosheets.  Note  that  a  static 
LLIP  method  (without  mechanical  agitation  or  ultra¬ 
sonication)  produces  larger  and  thicker  PCBM  crystals, 
similar  to  those  of  a  previous  report  (see  Figure  51).^ 
These  large  crystals  can  be  exfoliated  into  thinner 
nanosheets  by  ultrasonication  (Figure  SI),  but  by  this 
method,  the  nanosheet's  quality  (shape  and  size 
homogeneity)  is  compromised.  The  ultrathin  PCBM 
sheets  were  transformed  to  uniform  ID  PCBM  nano¬ 
structures  by  an  extended  ultrasonication  (30  min), 
which  results  in  a  folding  process  of  the  previously 
produced  PCBM  nanosheets  (see  Figures  1b  and  S2). 

The  resulting  material  contains  a  large  amount  of 
PCBM  rod-like  structures  (from  here  on,  denoted  as 


nanorods),  concomitant  with  some  small-sized  PCBM 
sheets,  following  the  sonication  (see  Figure  1c). 
Figures  2a  and  S3  reveal  that  the  nanorods  consist  of 
parallel  nanoribbons  in  different  numbers  with  clear 
boundaries  between  them.  Interestingly,  all  individual 
PCBM  nanoribbons  have  a  width  of  3.8  =b  0.3  nm 
(obtained  by  TEM),  which  exactly  matches  the  width 
of  four  aligned  PCBM  molecules.^^  Furthermore,  in  all 
nanoribbons,  the  molecules  along  the  nanoribbon 
edges  form  distinct  zigzag  patterns  (Figures  2b  and  S4), 
resembling  graphene  nanoribbons  but  in  molecular 
form.  The  lengths  of  the  nanorods,  on  the  other  hand, 
show  a  greater  distribution,  ranging  from  20  to  400  nm. 

Figure  2c  displays  an  individual  folded  nanoribbon 
(shown  by  arrows)  demonstrating  the  structural  flex¬ 
ibility.  In  contrast,  the  nanoribbon  in  Figure  2d  displays 
a  kink  at  two  locations  along  its  growth  direction.  This 
sharp  bend  indicates  the  presence  of  structural  de¬ 
fects,  and  in  this  case,  the  creation  of  two  separated 
divacancies  leads  to  a  formation  of  two  pentagonal 
rings  modifying  the  growth  direction  by  ~12°,  as 
illustrated  in  Figure  2d.  Another  important  observation 
from  HRTEM  analysis  is  that  adjacent  nanoribbons 
could  not  be  in  focus  at  the  same  time,  indicating  that 
they  are  located  at  different  heights  or  have  different 
orientation. 

We  note  that  during  the  synthesis  process  an  in¬ 
crease  in  the  sonication  time  results  in  improved  yields 
of  PCBM  nanorods;  likely  the  ultrasonication  process 
provides  necessary  energy  to  overcome  the  energy 
barrier  associated  with  the  folding  and  subsequent 
nanoribbon  formation.  However,  excessive  sonication 
power  or  time  breaks  the  PCBM  nanosheets  and  also 
fragments  the  already  formed  PCBM  nanorods.  Addi¬ 
tionally,  due  to  the  strong  hydrophobic  nature  of 
PCBM,  the  polarity  of  the  solvent,  used  during  the  LLIP 
and  ultrasonication  processes,  should  play  a  key  role  in 
the  formation  of  the  PCBM  nanorods.  To  examine  this 
effect,  the  polarity  of  the  reaction  environment  during 
the  synthesis  process  (LLIP  step)  was  decreased  by 
replacing  methanol  with  ethanol  (where  ethanol  ex¬ 
hibits  a  smaller  dielectric  constant).  The  resulting 
material  mostly  consists  of  PCBM  nanosheets,  but  after 
subjecting  the  sample  to  the  same  sonication  proce¬ 
dure  as  previously  described,  we  rarely  observed  any 
ID  nanostructures.  Interestingly,  adding  water  to  the 
system  (water/ethanol  =  2:3)  and  then  sonicating  the 
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Figure  2.  TEM  micrographs  of  ID  PCBM  nanostructures,  (a)  High-resolution  TEM  image  of  nanorods,  comprising  stacks  of 
individual  nanoribbons,  (b)  Larger  magnification  displays  that  each  nanoribbon  is  four  PCBM  molecules  wide  (3.8  ±  0.3  nm) 
packed  in  hexagonal  arrangement,  (c)  Folded  nanoribbon,  visualized  by  a  simulated  nanoribbon  at  the  right, 
(d)  Defective  nanoribbon  where  changes  in  its  growth  direction  (a  ~  12-14°)  are  caused  by  separated  two  divacancies, 
generating  two  pentagonal  member  rings  (marked  yellow  in  (d)).  The  dots  in  (b)  and  (d)  are  a  superposition  of  the  PCBM 
molecules  from  the  HRTEM  images. 


sample  once  again  led  to  the  formation  of  a  large 
amount  of  PCBM  nanorods,  suggesting  that  the  main 
driving  force  is  related  to  the  solvent— PCBM  interface. 
Note  that  the  yield  of  nanorod  formation  is  justified 
based  on  the  large  number  of  TEM  images  from  dif¬ 
ferent  batches  of  similar  samples.  Indeed,  adding  water 
to  samples  prepared  with  methanol  (water/methanol  = 
1:3)  significantly  increased  the  nanorods'  yield  further, 
and  thus  a  trend  regarding  the  nanorods'  yield  follows: 
ethanol  <  water/ethanol  <  methanol  <  water/methanol. 

Electron  energy  loss  spectroscopy  (EELS)  on  PCBM 
sheets  shows  no  signal  from  chlorine,  indicating  that 
the  PCBM  nanosheets  do  not  exhibit  solvent  molecules 
inside  their  structures.  The  EELS  data  are  supported 
by  thermogravimetric  analysis  (TGA)  data,  which  show 
only  a  weight  loss  of  less  than  1%  upon  heating  the 
PCBM  nanoribbons  to  300  °C  in  nitrogen  atmosphere 
(Figure  S5).  Additionally,  we  performed  a  thermal  an¬ 
nealing  on  PCBM  nanosheets  (200  °C  in  vacuum  for 
12  h),  and  the  obtained  material  did  not  exhibit  any 
structural  change,  such  as  cracks  or  cavities  that  are 
commonly  observed  in  solvated  fullerene  structures 
due  the  evaporation  of  trapped  solvent  molecules.^'^^ 
The  absence  of  structural  changes  upon  heating  to¬ 
gether  with  the  EELS  and  the  TGA  data  clearly  suggests 
that  the  structures  do  not  contain  any  solvent  mole¬ 
cules  and  that  the  PCBM  nanosheets  are  composed 
of  bilayers  of  two  single-molecular  layers  of  PCBM 
stacked  on  top  of  each  other.^  Such  a  system  also  has 
a  formation  energy  lower  than  that  of  other  possible 
configurations,  as  explained  below. 

Energetic  Stability  of  PCBM  Sheets  and  Nanoribbons.  We 
performed  density  functional  theory  (DFT)  calculations 


to  understand  the  molecular  structure  of  PCBM  nano¬ 
sheets  and  the  subsequent  nanoribbon  and  nanorod 
formation.  The  crystalline  structure  of  PCBM  has  not 
been  clearly  identified,  and  the  few  existing  reports 
suggest  that  the  PCBM  crystal  structure  depends  on 
the  growth  conditions.^^'^®'^^  Zheng  etal?"^  performed 
a  detailed  study  and  suggested  that  the  PCBM  crystal¬ 
lizes  in  the  hexagonal  lattice,  which  agrees  with  our 
HRTEM  observations  (see  Figure  2),  and  we  chose  such 
structures  as  a  starting  point.  Additionally,  we  investi¬ 
gated  cubic  systems  generated  from  a  simple  cubic 
crystal  structure  (see  Figures  56  and  57).  A  total  of  12 
hexagonal  and  square  monolayers  and  bilayers  were 
studied,  and  our  results  indicate  that  the  most  energe¬ 
tically  favorable  structure  comprises  a  hexagonal  PCBM 
bilayer  with  an  A— B  stacking  pattern  with  the  side  chains 
pointing  outward.  The  optimized  structure  of  the  hex¬ 
agonal  A— B  bilayer  (Ha,b-1  b)  and  the  binding  energy  (Eb) 
as  a  measure  of  energetic  stability  are  shown  in  Figure  3. 

The  existence  of  the  H-bonds  in  PCBM  crystals  has 
been  previously  observed^®'^^  and  plays  an  important 
role  in  the  stabilization  of  PCBM  crystals.  We  observe 
such  effects  on  the  studied  monolayers  and  bilayers, 
where  systems  with  a  H-bond  network  are  highly 
beneficial,  compared  to  configurations  where  the 
H-bonds  are  not  present  or  are  scarce  (see  Figure  3c). 
The  H-bond  network  is  especially  promoted  when  the 
side  chains  are  directed  outward  from  the  layer.  In 
fact,  the  H-bonds  and  the  improved  packing  strongly 
promote  the  formation  of  hexagonal  arrangement. 
As  previously  stated,  the  HA,B-1b  is  the  most  stable 
system,  suggesting  that  the  PCBM  molecules  first  form 
hexagonal  bilayers  with  an  A— B  stacking,  and  then  by 
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Figure  3.  (a,b)  Top  and  side  view  of  a  hexagonal  PCBM  bilayer  with  an  A-B  stacking  (labeled  Ha,b-1  b).  White  dots  represent 
hydrogen  atoms,  gray  dots  represent  carbon  atoms,  and  red  dots  represent  oxygen  atoms.  The  most  favorable  configuration 
exhibits  side  chains  pointing  outward,  which  helps  to  minimize  the  contact  of  Ceo  with  the  polar  solvent.  The  interlayer 
distance  is  ~8.1  A  as  determined  by  dispersion-corrected  DFT.  (c)  Energetic  stability  of  12  different  square  and  hexagonal 
PCBM  monolayers  and  bilayers.  Negative  values  indicate  energetically  favorable  structures.  The  separation  on  the  x-axis  is  for 
clarity  only.  Detailed  models  of  different  configurations  are  shown  in  Supporting  Information  (Figure  S7). 
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Figure  4.  Schematic  procedure  for  PCBM  nanorod  formation,  (a)  Sonication  during  the  liquid-liquid  interfacial  precipitation 
process  creates  bilayers  of  PCBM,  where  the  side  chains  are  protecting  the  Cgo  from  the  polar  solvent.  (b,c)  Further  sonication 
promotes  bilayered  PCBM  sheets  to  fold  themselves  while  removing  the  trapped  solvent  at  the  interfaces.  (d,e)  Continued 
sonication  cuts  the  bilayers  and  creates  smaller  micelle-like  units,  called  nanoribbons,  which  are  agglomerated  and  seen  as  a 
nanorods. 


ultrasonication,  the  folding  and  the  subsequent  nano¬ 
ribbon  and  nanorod  formation  are  performed.  Here,  it 
is  noteworthy  to  mention  that  the  difference  from  a 
bulk  PCBM  with  a  hexagonal  close-packed  (hep)  crystal 
structure  and  the  proposed  hexagonal  bilayer,  Ha,b-1  b, 
is  the  orientation  of  the  side  chains.  The  Ha,b-1  b  system 
exhibits  a  micelle-like  configuration,  where  the  more 
hydrophilic  tails  protect  the  hydrophobic  Ceo  the 
polar  solvent,  which  we  believe  in  our  case  is  a  key 
point  for  the  nanoribbon  formation.  In  this  way,  the 
stacking  of  bilayers  generates  a  new  intermolecular 
distance  of  ~15  A  (see  Figure  S8);  this  interbilayer 
distance  might  be  responsible  for  an  additional  peak 
observed  in  the  X-ray  diffraction  (XRD)  analysis  of 
PCBM  nanosheets.  The  XRD  pattern  shows  a  distinct 
peak  at  29  =  5.6,  while  the  rest  matches  the  pattern 
obtained  from  bulk  PCBM  powder  (see  Figure  58). 

The  nanoribbon  formation  is  probably  generated  by 
subdividing  the  bilayer  and  creating  smaller  micelle-like 


units,  being  four  PCBM  molecules  wide,  where  the  oxygen¬ 
ated  side  chains  protect  the  hydrophobic  Ceo-  The  nano¬ 
rods  are  then  formed  by  stacking  multiple  nanoribbons. 
A  simple  schematic  procedure  is  described  in  Figure  4. 

Now,  considering  the  HA,B-1b  bilayer  structure  ob¬ 
tained  by  DFT,  we  created  PCBM  nanoribbons  with 
zigzag  edges  similar  to  the  ones  observed  in  the 
HRTEM  in  Figures  2  and  54.  The  zigzag  PCBM  nano¬ 
ribbons  grow  along  the  (1 1  —20)  crystal  direction  with  a 
hep  unit  cell.  According  to  classic  molecular  dynamics 
(MD)  simulations,  these  PCBM  nanorods  are  stable  at 
room  temperature  (both  in  vacuum  or  in  a  methanol 
solution),  and  when  methanol  is  introduced,  a  micelle¬ 
like  configuration  is  formed.  The  PCBM  nanorod  is 
depicted  in  Figure  5  and  Figure  59.  Finally,  we  simu¬ 
lated  a  TEM  image,  using  the  multislice  approach  as 
described  by  Gomez-Rodriguez  et  from  the  top 
view  of  a  nanoribbon  (inset  in  Figure  5c),  with  a  very 
good  match  with  the  experimental  TEM  image. 
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Figure  5.  (a)  Top  view  and  (b)  cross  section  of  a  PCBM  nanoribbon  in  solution  obtained  from  molecular  dynamics.  Note  that 
in  (b),  the  side  chains  cover  the  PCBM  nanoribbon  in  all  directions.  Solvent  molecules  and  H  atoms  are  omitted  for  clarity, 
(c)  HRTEM  micrograph  of  a  PCBM  nanoribbon;  the  inset  shows  a  simulated  TEM  image  of  (a). 


Based  on  our  experimental  and  theoretical  observa¬ 
tions,  we  suggest  that  the  nanoribbon  formation  is 
promoted  by  the  strong  and  continuous  sonication 
process,  where  well-ordered  PCBM  bilayers  are  first 
damaged,  creating  fissures,  defects,  and  eventually 
breaking  them  into  smaller  sections.  The  sonication 
process  promotes  a  reconstruction  until  the  smallest 
micelle-like  units  are  obtained,  in  this  case,  being  just 
four  PCBM  molecules  in  width.  Therefore,  the  process 
selectivity  is  a  consequence  of  thermodynamic  limita¬ 
tions,  where,  beyond  this  limit,  the  system  becomes 
energetically  unfavorable  and  promotes  the  self-as¬ 
sembly  of  PCBM  molecules  into  nanoribbons.  A  simple 
schematic  procedure  is  described  in  Figure  4.  On  the 
other  hand,  wider  PCBM  nanoribbons  might  not  be 
observed  simply  because  the  ultrasonication  process 


continuously  breaks  them  into  smaller  units,  where 
most  of  them  reach  the  lower  limit  (four  PCBM 
molecules). 

CONCLUSION 

We  have  demonstrated  a  solution-based  method  for 
synthesizing  crystalline  ID  PCBM  nanostructures  with 
a  uniform  width  of  3.8  =b  0.3  nm,  equivalent  to  four 
PCBM  molecules,  and  lengths  of  20— 400  nm.The  nano¬ 
ribbons  show  well-defined  crystalline  structure,  com¬ 
prising  PCBM  molecules  in  a  hexagonal  arrangement 
without  trapped  solvent  molecules.  The  yield  of  the 
nanoribbons  increases  with  the  polarity  of  the  used 
solvent.  Given  the  dimensions  and  the  physical  proper¬ 
ties  of  the  PCBM  nanoribbons,  it  is  plausible  that  they 
should  be  ideal  as  acceptor  material  in  organic  solar  cells. 


MATERIALS  AND  METHODS 

Synthesis  of  PCBM  Nanoribbons.  Phenyl-Cei -butyric  acid  methyl 
ester  (Sigma-Aldrich;  99.9%)  was  dissolved  in  chloroform 
(concentration  of  1  mg/mL)  and  ultrasonicated  for  30  min  to 
prepare  a  homogeneous  dispersion.  The  liquid-liquid  interfa¬ 
cial  precipitation  method  was  used  to  produce  the  PCBM  nano¬ 
sheets.  Prior  to  sample  preparation,  the  PCBM-chloroform 
solution  was  heated  to  60  °C,  and  then  10  mL  of  methanol 
was  gently  added  to  1  mL  of  the  PCBM-chloroform  solution. 
The  mixture  was  then  directly  placed  in  an  ultrasonic  bath  and 
sonicated  for  5  min  (using  a  power  setting  of  9  on  a  USC300D 
ultrasonic  bath  from  VWR,  frequency  45  kHz  and  effective 
power  80  W).  During  sonication,  the  color  of  the  blend  changed 
from  fuzzy  brown  to  clear  brownish.  The  new  mixed  suspension 
was  then  stored  at  room  temperature  for  24  h  in  a  tightly  closed 
bottle.  The  PCBM  nanorods  and  nanoribbons  were  then  pro¬ 
duced  by  collocating  the  previous  sample  in  an  ultrasonic  bath 
for  30  min  (same  power  setting). 

Sample  Characterization.  The  prepared  samples  were  trans¬ 
ferred  on  TEM  grid  by  dipping  the  grid  in  the  dispersion 
containing  PCBM  structures.  Note  that  before  TEM  grid  pre¬ 
paration  the  sample  was  collected  by  centrifugation  and  redis¬ 
persed  in  fresh  methanol.  The  TEM  measurements  were  done 
using  a  JEOL  2010  (operated  at  80  keV  with  a  LaBs  gun),  JEOL 
2100  (operated  at  200  keV),  or  JEOL  1230  (operated  at  80  keV). 
For  TEM  characterization  of  the  heat-treated  structure,  the 
samples  were  annealed  on  the  TEM  grid  at  200  °C  in  vacuum 
for  1 2  h. 

Ab  Initio  Computational  Details.  Ab  initio  computations  were 
performed  using  DFT^"^  and  the  generalized  gradient  approx¬ 
imation  with  the  Perdew,  Burke,  and  Ernzerhof  modeP^  as  the 
exchange-correlation  term.  The  wave  functions  for  the  valence 
electrons  were  represented  by  a  linear  combination  of  pseudo- 
atomic  numerical  orbitals^^  using  a  double-^  basis. The  real-space 


grid  used  for  charge  and  potential  integration  was  equivalent 
to  a  plane  wave  cutoff  energy  of  200  Ry.  The  integration  of 
the  Brillouin  zone  was  carried  out  using  an  8  x  8  x  8 
Monkhorst-Pack^^  grid  for  bulk  systems;  however,  a  1  x  5  x 
5  Monkhorst— Pack  grid  was  used  for  2D  systems.  A 
Methfessel- Paxton^®  smearing  of  25  meV  was  employed  to 
aid  convergence  on  the  self-consistent  electron  density.  Peri¬ 
odic  boundary  conditions  were  used  with  at  least  30  A  of 
vacuum  to  avoid  intercell  interactions.  All  systems  were  fully 
relaxed  until  the  maximum  force  was  <0.04  eV/A.  The  interlayer 
distance  was  calculated  considering  Grimme's  dispersion 
corrections.^^  The  DFT  calculations  were  performed  using  the 
SIESTA  code.""® 

We  studied  diverse  square  and  hexagonal  monolayers  and 
bilayers.  The  square  systems  were  built  from  a  bulk  PCBM  with  a 
simple  cubic  unit  cell,  a  =  9.8  A.  The  hexagonal  systems  were 
obtained  from  a  bulk  PCBM  with  a  hep  crystal  structure,  where 
a  =  b  =  9.8  k  and  c  =  23.  A.  Finally,  the  energetic  stability  of  at 
least  1 2  different  structures  were  compared  by  calculating  their 
binding  energy.  The  binding  energy  was  defined  as  Eb  =  ^pcbm  - 
^system/  where  EpcBM  is  the  energy  of  the  isolated  PCBM  molecule 
and  Esystem  is  the  energy  of  the  bulk,  bilayer,  or  monolayer 
system,  accordingly. 

Molecular  Dynamics  Computational  Details.  The  MD  simulations 
were  carried  out  using  the  large-scale  atomic/molecular  mas¬ 
sively  parallel  simulator"^^  and  ReaxFF"^^'^^  as  a  reactive  force 
field.  For  the  ReaxFF  simulations,  we  used  the  force  field 
parameters  reported  by  Rahaman  et  these  parameters 
have  been  shown  to  correctly  represent  organic  molecules  in 
solvents  such  as  water  and  alcohols.  MD  was  carried  out  under  a 
constant  number  of  atoms  and  volume,  while  the  temperature 
was  controlled  by  a  Berendsen  thermostat^^  with  a  0.1  ps 
damping  constant;  these  conditions  are  denoted  as  NVT.  First, 
an  energy  minimization  was  performed,  and  then  the  tempera¬ 
ture  was  increased  by  a  constant  rate  equal  to  5.0  K/ps  up  to 
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298  K;  once  the  temperature  was  reached,  the  MD  simulation 
was  continued  for  another  200  ps  using  a  time  step  of  0.25  fs. 
The  nanoribbons  were  based  on  the  hexagonal  bilayer  with  an 
A— B  stacking  and  the  side  chain  pointing  outward.  Its  ID  unit 
cell  contained  16  PCBM  molecules;  however,  the  studied  super¬ 
cell  consisted  of  48  PCBM  molecules  (4224  atoms),  where  the 
periodic  dimension  had  a  lattice  parameter  of  34.64  A,  while 
the  other  two  were  kept  at  66  A  to  avoid  lateral  interactions. 
The  MD  simulations  of  nanoribbons  in  methanol  solution  were 
performed  using  similar  conditions,  but  now  randomly  packing 
2260  methanol  molecules  around  the  PCBM  nanoribbons 
(1 8  776  atoms  in  total). 
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